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Abstract 



We study the pion electromagnetic and 7* + tt° — ► 7 transition form fac- 
^ ■ tors at intermediate momentum transfer. We calculate soft, nonperturbative 

corrections to the leading perturbative amplitudes which arise from the qq- 
component of the pion wave function. We work in Minkowski space and use a 
Lorentz covariant, gauge-invariant generalized perturbative integral represen- 
tation for the qq amplitudes. For the transition form factor we find relative 
insensitivity to the detailed nonperturbative structure of the wavefunction for 
\q 2 \ > 10 GeV 2 , whereas considerable sensitivity is found for the electromag- 
netic form factor. 
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I. INTRODUCTION 



There has recently been considerable attention focused on the applicability (or other- 
wise) of the leading twist perturbative QCD (pQCD) analysis in descriptions of exclusive 
processes at large momentum transfers. Theoretical studies of such processes are typically 
justified on the basis of the factorization of the hard, perturbative and the soft, nonper- 
turbative contributions to hadronic matrix elements. This factorization appears explicitly 
through the operator product expansion (OPE) M, where amplitudes are determined by 
the diagonal or almost diagonal matrix elements of nonlocal operators. It is applied to 
far-off-diagonal (i.e., large momentum transfer) matrix elements of local operators in a par- 
ton model description of exclusive hadronic form factors [|2],|3|,|4]|| . The largest discrepancy 
between experimental data and the leading order perturbative calculations is observed in 
the predictions for absolute normalization of hadronic contributions to cross sections such 
as the K factor in the Drell-Yan process and electromagnetic form factors. This may be 
due to large corrections from next to leading order (in a s ) terms and/or contributions from 
operators with higher dimension (twist). For example, the pion electromagnetic (e.m.) form 
factor [F(Q 2 )] at large spacelike momentum transfer has been the subject of several studies, 
which showed that there were difficulties associated with the application of leading order 
pQCD arising from the appearance of intermediate states with small virtuality. These arise 
from the end-point region in the loop integration over the quark momenta [||],|6],|7| in what 
are necessarily assumed to be hard amplitudes. In particular, they may lead to large pertur- 
bative corrections, rendering the perturbative expansion inconsistent. Recent analyses |||| 
have shown that the inclusion of Sudakov-type corrections and transverse momenta in the 
analysis can suppress the unwanted end-point contributions and lead to a self consistent 
perturbative expansion for momentum transfers as low as a few GeV. There is, however, an 
open question on the role of the higher twist operators. This point was first emphasized 
in Ref. |J, and to the best of our knowledge for the case of the pion e.m. form factor 



has only been examined in a model independent way in Refs [plj and ||11||. In Ref. (10 
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the contributions from twist-3 operators in a parton model factorization approach has been 
calculated and shown to be even larger than the naive asymptotic one for Q 2 < 20 GeV 2 . In 
Ref ]TI] the behaviour of the form factor at Q 2 ~ 10 GeV 2 has been calculated within the 
QCD sum rule approach |T!| and was also shown to be dominated by contributions from 
higher dimension operators. In this work we wish to address the question of the role of 
nonperturbative corrections to descriptions based both on parton model factorization and 
the OPE. In Sec. we use a Lorentz covariant wavefunction formalism to extend the parton 
model approach to the pion electromagnetic form factor. This is an extension of initial work 
presented elsewhere |13|- In Sec. [HI] we calculate the nonperturbative corrections to the 
7*7r — > 7 form factor [T(Q 2 )} as an example of an application of the OPE techniques. We 
give our conclusions in Sec. |TV]. 



II. PION ELECTROMAGNETIC FORM FACTOR 

A complete representation of the pion e.m. form factor, F(Q 2 ) is given in Fig. [I] and 
involves a pointlike coupling of the photon to a quark and anti-quark, (all dressings etc. are 
contained in M M ); M M = F(Q 2 )(P^ + P >fl ), represents the matrix element of the e.m. quark 
current, J^(0) = : 4 , (^) e ql' J 'i , (^) between pion states of momenta P and P + q = P', 
(Q 2 = — q 2 ). In QCD, as the momentum transfer, Q 2 — ► oo, the leading contribution to the 
form factor factorizes into a convolution of a hard scattering kernel (T) which involves a 
minimal number of hard, highly virtual parton lines, with soft distribution amplitudes (0) 
for the initial and final pion. Assuming factorization, the qq state contribution to the form 
factor in M M can be written as PH], 

F(Q 2 ) - F(Q 2 ) pQCD = £ dxdy^ A (x, \ 2 )T(x, y; Q 2 , \ 2 )Mv, A 2 ) , (1) 

where A refers to the factorization scale and where the subscript A on 0^ indicates the 
soft pion amplitude arising from an axial-vector spinor structure. In pQCD, large and 
small virtuality contributions from loops are factorized into T and respectively, x and 
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y are fractions of momentum relative to the momentum of the initial and final on-shell 
massless pion respectively. The hard scattering amplitude to leading order in a s is given 
by a Born diagram with a one gluon exchange projected on a collinear, massless, on-shell 
qq state with an axial- vector spinor structure in initial and final channels (i.e., the twist- 2 
contribution). Similarly, the corresponding distribution amplitudes <pA are obtained by a 
collinear projection of the axial-vector ($a) component of the matrix element, 



In Eq. (0) sums over the color indices i and j are understood and the parallel transport 
operator can be written as 



where V denotes the path-ordering operator and a and (3 are spinor indices. The color 
labels i and j in Eq. (0) have the ordering shown since the parallel-transport operator usually 
appears in the combination ip(y)P(y, x)ip(x), where the matrix summation over color indices 
is implied. The gluon field operators are denoted by A£(x), where a = 1, 2, . . . , 8 and where 
A^(x) = (l/2)\ a A%(x). g s is the strong coupling constant. In the amplitude $ the incoming 
pion has 4-momentum P, the quark has 4- momentum p + (P/2), and the antiquark has 
4- momentum p — {P/2). We see that the parallel transport operator plays the same role 
in $ as that played by the gluon flux tube in the valence quark model in that it ensures 
local color gauge invariance fl4| . The color-singlet amplitude $ is identical to the familiar 
Bethe-Salpeter (BS) amplitude when the q and q couple at a single point, (i.e., when z = 0). 
This is just the situation that occurs, for example, in a calculation of the pion decay constant 
/„-, where there is a coupling to a pointlike W or Z, (see later). 

We shall estimate the size of 0(1/ Q 2 ) corrections to the asymptotic behaviour -F p q CD of 
the form factor once the collinear projection has been removed; these corrections are due to 
the presence of nonzero quark virtualities in the soft &a part of the wave function in Eq. (0) , 






(3) 



where 
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P) = 75 $p(p, P) + P 75 $ A (p, P) +#75^(P, P) + fe/,P] 7 5$T(p, ^) • (4) 

Here $j(p, P) for i = P, A, A', and T are scalar functions of p 2 , (p • P), and P 2 . The 
subscripts P, A, A', and T denote the spinor matrix structure of the particular component 
of $. The three functions $p, $>a, and $t are even functions of p ■ P, whereas $,4/ must be 
an odd function of p ■ P. We will work here in the chiral limit where the pion is massless 
and so for an on-shell pion we have P 2 = M 2 = 0. 

A standard assumption [15j is that the scalar functions $i(p, P) can be written in terms 



of a generalized perturbation-theory integral representation, which for equal mass particles 
we can express as 



/l roc 
dx / dv 2 
-1 Jo 



00 ! 

E 

/3=2 ' 

7; 1 1?' (5) 

(p-xiPJ - [z/ 2 + m 2 - |(1 - x 2 )P 2 J +ze 

where m is the u and quark mass (assumed degenerate for simplicity) and the g are scalar 
functions generally taken to be polynomially dependent on p M ||15|| . In the above (3 is an 
integer > 2. Once the perturbative corrections are taken into account, the renormalized 
parameters and the scalar functions g pick up a log A 2 dependence, where the factorization 
scale A divides soft and hard contributions. We shall return to this point later. In order to 
ensure confinement in such a representation we can require, for example, that the functions 
g have no support [i.e., g(x,u,p, P) = 0] in the region where v 2 < (P 2 /4) — m 2 . In the 
chiral limit for on-shell pions (m = 0) the surviving nonperturbative scale is that typical of 
dynamical chiral symmetry breaking and confinement. We will denote this scale here by fi 
and we clearly expect fi < 1 GeV, at A ~ 1 GeV. 

With this motivation in mind we now make an assumption about the form of the functions 
g~p, i.e., we assume 

gpi(x,v;p,P) = 5(v 2 - /j 2 ) N pi g pi (x) , (6) 
where the mass scale fi, and the functions gpi{x) have yet to be specified. The functions 



will be chosen to be dimensionless and so the normalizations Npi are included for dimen- 
sional reasons, and are to be determined by fitting the pion wave function to some physical 
quantities, (e.g., f n as discussed later). With this choice we then find from Eq. (||) that 

As already pointed out, only the soft part of the amplitudes $j(p, P) are to be used in the 
calculation for the form factors and has to contribute only at small virtualities, (p± \P) 2 < 
A 2 , with the contributions from hard virtualities being factorized into the hard scattering. In 
order for that to happen, we require that the soft part of the amplitude $j given by Eq. (|7|) 
has its short distance behaviour softer than the one corresponding to a high momentum 
gluon exchange which leads to the hard scattering amplitude. In other words, any large 
relative momentum behaviour in Eq. (^) that would be identical to that from a hard gluon 
exchange has to be subtracted, in order to avoid double counting. 

In order to isolate that part of P) which corresponds to the general light-cone 
z 2 < 1/A 2 short distance behaviour of the matrix element in Eq. (0), one has to study 
the p\ > A 2 behaviour of the corresponding light cone projection <3> L (p + ,p±, P) of &a, 
since the dominant asymptotic behavior arises from this. We use standard notation, where 
z ± = z° ± z 3 , zj_ = (z 1 ,^ 2 ), etc., and work in the infinite momentum frame where P — > 
(P + , 0~, 0j_). As usual, we also define the fraction of relative longitudinal momentum as 
x = 2p + /P + , [we will see shortly that this x can be related to the x in Eq. (0)]. The light 
cone wave function is defined as 

$ L( V, P± ,P) = $a{p,P) ■ (8) 



Substituting the expression given by Eq. (0) for &a in Eq. (|8|) we obtain for p\ > A 2 , 

(9) 



$ LC V,p ± ,P) oc V N ^9pa(x) _N 2A g 2A (x) 



5(3,2 + O 



0=1 p+ (pi + fi 2 )"- 1 ~ p\p + r" 1 ~ \ x2 pi, 

On the other hand the pQCD expansion of $ LC for p\ > A 2 corresponding to a hard gluon 
exchange is given by JL6J 
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( n+1 1 2 \ 

7 n = 2C F 1 + 4^-- ,7f = CV, (10) 

V (n + l)(ri + 2)y 

where Cp = 4/3 and /?o = 11 — |n/ with n/ being the number of quark flavors. The coeffi- 
cients a£> contain all of the nonperturbative information and are sensitive to the behaviour 
of the soft LC wave function i.e. for p\ < X 2 . The functions G^/ 2 (x) are the usual Geigen- 
bauer polynomials of order 3/2 and contain all of the longitudinal momentum information. 

An examination of the double integral (over x and p~) encountered in deriving Eq. @ 
shows that the ^"-integral will give zero unless we have x = 2p + / P + . Hence as previously 
stated we can identify the x appearing in the functions gpi{x) with the fraction of longitudinal 
momenta. Comparing the leading terms in the \ 2 jp\ expansion of Eqs. (|^) and ([H]), we 
see that leading power behaviour of the (3 = 2 term in Eq. (^|) is the same as the one given 
by a hard gluon exchange tail in Eq. ([II]). Thus the (3 = 2 term should be excluded from 
the soft pion wave function in this limit, as explained above. Furthermore, comparing the 
x-dependence of the corresponding asymptotic terms, we see that if logarithmic corrections 
of Eq. fllPf ) were taken into account in Eqs. (0) and (H), one would have for p\ > A 2 , 

92 A (x) ^ g 2A (x,pl) oc (1 -x 2 )E4 A) f !°^M? CD M GT{-) , (11) 

n \ l0 gl A _L/ A QCD) ) 

On substituting Eq. (|7D into Eq. fllPP we also see that the coefficients 's can be expressed 
in turn in terms of the (normalized) Geigenbauer moments of the functions gpA{x), for 
which gpA^x) —> gpA{.x,P±) for p\ < A 2 . This relationship can then be inverted using the 
orthonormality of the Geigenbauer polynomials [with the measure (1 — x 2 )} to give the g@A 
in terms of the coefficients . The (3 = 2 component in the evaluation of the coefficients 
(for p\ < A 2 ) must be excluded in Eq. ([U]) because it would generate a nonperturbative 
model-dependent cutoff, [i.e., a factorization scale (A) dependence], for physical amplitudes 
like f n for example. Thus, in keeping with the usual concept of factorization we have for 



(3 > 3 

9Pa(x) -> gp A (x, A 2 ) oc (1 - x 2 ) £ a^GT{x). (12) 

n 

Hence, from Eqs. (0) and ( |TTD it follows that g 2 A(x, p 2 ± ) should be identified with the pQCD 
evolved quark distribution amplitude. For the nonperturbative (/3 > 3) components of the 
wave function (/3 > 3) we will assume here that such terms can be adequately represented 
by the piece that will dominate the X 2 /p 2 ± expansion, i.e., we neglect (3 > 4 and use a (3 = 3 
piece alone. This is not essential, but obviously simplifies the analysis without altering the 
asymptotic behaviour. We shall comment further on this point later. 

The functions <?3a(^, A 2 ) can then be assumed to correspond to typical soft quark distri- 
bution amplitudes. Following the above arguments we will simplify the notation. We shall 
now understand that $;(p, P) of Eq. (^) refer to only the soft nonleading contributions with 
(3 = 3. Hence we write in place of Eq. (|7D for $j 

*Mm-»»£ t ^_$*'tl + ur (13) 

The distribution amplitudes are typically normalized such that 

dx gpi(x, A 2 ) = 1, (14) 



and in particular for p\/\ 2 — > oo the usual asymptotic result for g 2 A is 

3 

g 2 A(x) ^ g 2A (x,oo) = -(I- x 2 ) . (15) 

while the soft quark distribution amplitude might, for example, be taken to be parameterized 
in terms of the Chernyak and Zhitnitsky moments 

g 3A (x, A 2 ) -> g 3A (x, ~ 1 GeV 2 ) ~ g cz (x) = ^(1 - x 2 )x 2 . (16) 

The normalization constant iV^ can be fixed by the fact that the 0-th moment of $ is 
directly related to the pion decay constant f n = 93 MeV. Only $^ and can contribute 
to fn and we shall for convenience make a further simplifying assumption, which is that 

8 



nonperturbative piece can be neglected. A straightforward calculation of f n then only 
involves $ A and gives 

U = -A- t dx[N 3A g 3A (x, A 2 )] . (17) 
on n J-i 

It follows from Eq. (|T3]) that N 3A = — (87r 2 /3)/ 7r/ u 2 . This leads to our final expression for 
^ / m 8tt 2 / w r 1 , g 3A (x,\ 2 ){i 2 

®a(p,P)(\) = — I dx— ^ . (18) 

p — x\P) — /i 2 + ie 



It is worth noting that Eqs. (§), (|10D , and (|14"D guarantee that /„- is factorization scale (i.e., 
A) independent even when the pQCD corrections are included. If we were to generalize to 
a wave function containing (3 > 3 terms, then the unknown structure functions, g@i(x, A 2 ), 
(/3 = 4, • ■ ■) and corresponding scales, fipi, would need to be fixed from the normalization of 
additional physical matrix elements. Matrix elements involving higher number of covariant 
derivatives such as 

<(#(0)(i S) VtbC* 5r • • • (i B)^V(0)|P)(A) = ifAkl) A (x n ) G P» ■ ■ -P^ (19) 
involve the moments of the distribution amplitude gc(x, A 2 ) ||, 

^(x,A 2 )P^ = J ^•P)e-^^0|^(^)^(0) 7 ^(-^)| J P)(A), 

(x") G = / dxx n g G (x,\ 2 ) . (20) 



These are associated with the quark longitudinal momentum distribution in the presence of 
an additional gluon operator G^ u = \e^ ua pG a p. As described in the introduction, however, 
we restrict our analysis to the contribution from the qq sector only and thus are ignoring 
possible higher Fock sector distribution amplitudes. 

Returning now to the e.m. form factor calculation, the standard quark counting rules and 
the pQCD formula of Eq. <g) to leading order in Q 2 predict F(Q 2 ) -»■ F pQCB {Q 2 ) ~ 1/Q 2 
as Q 2 — > oo. The 1/Q 2 behavior comes from diagrams with a single gluon exchange in a 
quark loop. In order to maintain QCD gauge invariance, once gluon degrees of freedom 
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have been explicitly introduced, one should include the soft spectator gluon diagram of 
Fig. |2|b) along with the hard gluon exchange diagrams of Fig. |2]a). However, only Fig. |2|a) 
contributes to the leading 1/Q 2 behaviour and so the gauge dependence enters as an 0(1/ Q 2 ) 
correction to the leading behaviour. In order to show that the formalism presented here 
gives the right asymptotic behaviour we calculate the diagrams of Fig. ^|a) using the $^4 
contribution and an arbitrary covariant gauge where the gluon propagator has the form 
D a /3 = [—ga/3 + (1 — a)k a kp/{k 2 + ie)]/(k 2 + ie). We shall refer to this contribution to the 
form factor as F A - This gives then (q 2 = —Q 2 < and P' = P + q) 

dH d 4 k 



i(p + p'YF A (Q 2 ) = -6 J 



(2tt) 4 (2tt) 



tr 



2 

P' 



+ (-±)*a(1 - P)T"(k, I, -P, -P')<S> A (k - y , -P') S- 1 ( ( ) 



(21) 



3 y ^ 2 

where (2/3) and (—1/3) are the u and d quark charges respectively and where for the pion 



17,18 



$(p, P) = T$ T (p, -P)^ 1 = $(p, P) . (22) 

The transpose acts on the spinor indices and T = i^y 1 ^ 3 = = T _1 is the time-reversal 
operator in the usual Dirac representation. The hard scattering amplitude T M describing 
the piece of the diagrams which do not involve the wave functions is given by, 

!*»(*, I, P, P') = g 2 s [i a S F (k + P) 7 1 ® [i p }D aP (k - I) 

+ [rfS F {l + P'h a ] ® [l P ]Dap{k - I). (23) 



With Sp being the perturbative quark propagator. Using the form for $^ given by Eq. (|18| 
one finally obtains, 

Fa(Q 2 ) = ^- j^dxg 3A {x,X 2 ) j\a/ 93A {a/,\ 2 ) 



Jo 



^ 2 Q 2 e(i-x) 



+ (x — > X 



e(i - x)f + f\ [e(i - oa - x)(i - *of + ^ 2 
1 + 



(24) 
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Once all of the one loop corrections to the Born diagram are taken into account and renor- 
malization at fi R ~ Q 2 is performed [in order to cancel potentially large logs of the form 
log(Q 2 ///|j)], the unrenormalized coupling, g 2 gets replaced by ATia s {Q 2 ). Similarly, logs 
of the form log(Q 2 /A 2 ), with A introduced in a loop calculation in order to factorize soft 
and hard contributions, are responsible for the evolution of the initial, soft distribution 
amplitudes g 3A (x,X 2 ) in a form that is given by Eq. (|TTD ||. Thus the expression for the 
renormalized form factor F A can be obtained from Eq. fl24|) by the following replacement, 

g 2 s g 3A (x, X 2 )g3 A (x, A 2 ) -> Ana s (Q 2 )g 2A (x, Q 2 )g 2A (x f , Q 2 ) [1 + 0{a s )\ . (25) 

The remaining, finite perturbative corrections in Eq. (p5| ) are actually large 0, however we 
are concerned only with the role of the nonperturbative corrections and so in the following 
we do not discuss the perturbative ones. The 0(fi 2 /Q 2 ) terms on the RHS come from the 
gauge fixing terms proportional to (1 — a) and we can neglect these contributions here. In the 
asymptotic limit the leading contribution from the integral in Eq. ( pip to the renormalized 
form factor is of the form F A {Q 2 ) — > F p q CD (Q 2 ) , where 



2 _ 64vra s (Q 2 )/ 2 

pQCDlV ) — - 



1 dx 92AM) 
-1 1 — X 



(26) 



9Q 2 

When we also substitute the asymptotic form for g 2A from Eq. (|I5|) we recover the well known 
leading- twist pQCD result @,@,||]. As mentioned in the Introduction the twist-3 operators 
were shown to give also a sizable contribution to the e.m. form factor at intermediate 
momentum transfers. These corrections correspond to the contributions from the $p and 



&t pieces of the our wavefunction and since they mix under renormalization []T0[ both have 



to be included. We use the following generalization of the collinear twist-3 wave function 



introduced in Ref. [ 1 1 that contains both P and T components, 



-1 [(p — x\P) 2 — jl 2 



We write this as a (3 = 3 piece because the extra power of momentum squared in the 
numerator cancels a corresponding power in the denominator. We use f_ x dxg 3 p(x) = 1 and 
the normalization coefficient N 3 p calculated from the matrix element 
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(0|^(0) 75 ^(0)|P) (A) = *v^^. 



(28) 



where (qq) x^iGeV ~ ~~ (250MeV) 3 is the quark condensate. The contribution F P {Q 2 ) to the 
form factor coming from the twist-3 wavefunction above can be obtained in a way analogous 
to the Fa{Q 2 ) one and is given by, 
8rra s (Q 2 )(qq) 2 {Q) 



F P (Q< 



9/ 2 



x 



1 — x , , 1 — x'\ 



1 A rl rl 

dxg2p(x,Q 2 ) / dx'g 2 p(x',Q 2 )S(x,x')x d£ dr/rj 
-1 J-l Jo Jo 

+ (x <-> a;') 



(29) 



with pura 



(qq) 



(Q) 



'a s (Q 2 )Y m) 



(30) 



and g2p(x,Q 2 ) being the twist-3 distribution amplitude obtained from the g 3 p(x,X 2 ) one 
through the pQCD evolution equation coming from the renormalization of the vacuum-to- 



pion matrix element of the twist-3 operator [10] 



-7n/2/3 



( n+l l 

7 n = 2C F ll + 4^- 

g 3 p(x,\ 2 )=E bi n ) Pn(x), 



(31) 



and P n being the Legendre polynomials. Finally 



a 



S(xx') = ( — 



C tp , ot s (a ) 



O l = XX'Q 2 + yU 



(32) 



is the Sudakov form factor |j^riTJ] . Although we have so far been neglecting the 0(a s ) 
perturbative corrections, the Sudakov form factor, S(a 2 /Q 2 ) is needed in Eq. ( |29"| ) to keep 
Q 4 Fp(Q 2 ) finite in the Q 2 — > oo limit as required by perturbation theory [|T0[ |. 

In our treatment the pion e.m. form factor is F(Q 2 ) = Fa(Q 2 ) + Fp(Q 2 ), where Fa is 
given by by Eqs. ([24]) and (^5|) and Fp is given by Eq. (^9|) . We have used a running coupling 
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a s = 4:7i/j3 \og(Q 2 / Aqcd) with Aqcd = 150MeV. For the pseudoscalar distribution ampli- 
tude an asymptotic form g2p(x, Q 2 ) — > gp(x, oo) = 1/2 has been used. For g 2 A we have used 
the Chernyak and Zhitnitsky moments of Eq. fllED, i.e., g2A(x,Q 2 ) ~ ^^(A 2 ) ~ gcz(x) on 
the basis of Eq. fl25|) and since the Q 2 dependence due to the pQCD evolution is small. This 
is certainly reasonable since the finite 0(a s ) corrections to the Born amplitude have been 
neglected as well @. In our calculation the parameter /z has been set to be fi — 338MeV 
which is a mass-scale typical of both confinement and dynamical chiral symmetry breaking 
and it has been chosen such that our pion wavefunction give a correct normalization of 
the 7*7r — > 7 form factor which we analyze in the next section. In Fig. |3| we have shown 
four curves, where the lower solid curve is the prediction resulting from the above argu- 
ments. For comparison we also show the result for -F p qcd defined in Eq. fl2"6|) when we use 
g2A{x,Q 2 ) = gcz(x), which is the upper solid curve. The dashed curves correspond to the 
solid curves, except that we have used the asymptotic distribution g2A(x,Q 2 ) — > g2A(x, oo) 
given in Eq. (|T5| ) rather than gcz- The dashed curve which increases at low Q 2 is the one 
which contains Fp. At very large Q 2 all results converge as is of course to be expected. With 
the asymptotic form for g 2 A neither the pQCD formula nor the form factor calculated here 
appear to reproduce the data. In the Chernyak and Zhitnitsky case, which is more realistic 
for a distribution amplitude at A ~ lGeV, the introduction of the soft corrections in the 
wavefunction significantly reduce the size of the form factor below the pure pQCD result in 
the intermediate regime (~ 5 to 20 GeV 2 ). In Fig. |] we illustrate the overall Q 2 dependence 
of the 0(1/Q 4 ) (at large Q 2 ) corrections to the pQCD form factor F p q C d- We attempt to 
do this by plotting the dimensionless combination, 

SF(Q 2 ) = FpQCD J g2) :Zt m + Fp(Q 2 ) ■ (33) 

The net effect is a substantial correction at Q 2 ~ lOGeV 2 for both the Chernyak and 
Zhitnitsky distribution amplitude (solid line) and the asymptotic one (dashed line). 



13 



III. 7*tt u -> 7 FORM FACTOR 

The full representation of the 7*7r — > 7 form factor amplitude T^Q 2 ) where Q 2 = — <? 2 
and g 2 is the mass of the virtual photon is shown in Fig. [|. The relevant hadronic piece of 



the corresponding amplitude is given by [20 



= ij d 4 ze-^ z ($\TJ,{z)JM\P) = ^u P *q p P a T(Q 2 ) , (34) 

For \Q 2 \ -^ooan operator product expansion (OPE) of the two currents in terms of local, 
gauge invariant operators can be used to calculate F{Q 2 ). We wish to extrapolate the calcu- 
lation of the form factor from the asymptotic region, where only the twist- 2 piece of the OPE 
of the product J^{z)J v {Q) and the light-cone wave function of Eq. ((T0|) contribute [pOp , down 
to the intermediate region of momentum transfer. Furthermore, we wish to do this while 
maintaining the valence (qq) contribution only. This requirement corresponds to keeping par- 
ticular contributions to the matrix element (0\TJ fJi (z)J u (0)\P). This includes contributions 
of the form S (d^ . . . d /lm ip'y' lm+1 'y!iD IIm+2 • • • Dum+n^ ~ traces) , (which are twist-2 pieces that 
lead to the known perturbative result ||20||), contributions from the twist-2 "trace" pieces 
that correspond to replacing one or more pairs of covariant derivatives D^D V by g^D 2 , and 
finally higher twist contributions obtained by contracting pairs of covariant derivatives with- 
out the spin-2 g*" tensors. For example, S (d^ . . . d^ip^+'^D^^ . . . D 2 . . . D^ m+n il?j is 
a twist-4 contribution. S represents the usual symmetrization of the space-time indices. 
Clearly, there are no explicit gluons or sea-quarks in these operators. The sum of all such 
operators is equivalent to the replacement |2l], 



(0|TJ M ( 2 )J,(0)|P) -> 2e AB , ffp [0"A(*)]T!r [ 7 p 7 5 (0|T^(z)P(^, 0)^(0)|P)W] (35) 

where P(z, 0) is the parallel transport operator and A(z) = l/(47r 2 ,2 2 — is) originates from 
the vector part of a free massless Dirac propagator. A simple description of this is that 
we have a purely valence pion wavefunction with a photon attached to each of the q and q 
legs and with a massless quark line connecting the two photon vertices. This is shown in 
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Fig. ^]. The replacement represented by Eq. ( |35|) corresponds to assuming that all nonleading 
contributions can be subsumed into the valence pion wavefunction, which we immediately 
identify with the $,4 component of the covariant wavefunction from Sec.[TT[ Thus we shall 
drop the index A on the wavefunction. The leading (hard) term dominates asymptotically 
and as in the case of the pion e.m. form factor [see Eq. (|I])] is obtained through a collinear 
projection of the wave function in Eq. (|3~5"1). From Eqs. (|T0|) , (|TTD , Cj34|) , and (|35|) the pQCD 
behaviour is given by T(Q 2 ) — > T p qcd(Q 2 ), where 

T pQ c»(Q ) = ^ l^^^r ■ ( 36 ) 
In general, for our soft wave function of Eq. (|18"D, we find that Eqs. and ( |35|) lead to 



3Q 2 J -1 1-x 



| 2/x 2 log (l + gr(l-x)-ie 



(37) 



Q 2 (l-x) 

which is the form that we use below the factorization scale (\Q 2 \ < A 2 ) with g 3 some soft 
quark distibution amplitude. Above the factorization scale (\Q 2 \ > A 2 ) the perturbative 
corrections to the c-number coeficient [A (z)} in Eq. (^) lead to the replacement 

2 3 0r,A 2 )^ 2 0r,Q 2 ) (38) 

through Eq. (|TTD with the coefficients given by Eq. (]12D . Thus, above the factorization 
scale we simply use the same distribution amplitude as below, but include the logarithmic 
correction factor of Eq. ([TT|) . 

In Figs. 0a) and b) we plot our predictions for the T(Q 2 ) form factor in the spacelike 
region (Q 2 > 0) and timelike region (Q 2 < 0) respectively. In the spacelike (timelike) re- 
gions the upper (lower) solid line is the pQCD prediction of Eq. (^) with the Chernyak 
and Zhitnitsky distribution amplitude used as an input for the soft distribution amplitude, 
i.e., (g 2 (x,Q 2 ) — > gcz{x)) for \Q 2 \ < A 2 = lGeV 2 . For \Q 2 \ > A 2 pQCD as explained above 
the evolution given by Eq. ( |TTD was used to generate perturbative logarithmic corrections to 
T(Q 2 ). The upper (lower) dashed curves represent the pQCD predictions in the spacelike 
(timelike) regions for the asymptotic distribution amplitude [see Eq. (|15|) 1 (for which pQCD 
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evolution has been neglected). Our prediction for the form factor [see Eq. using the 

Chernyak and Zhitnitsky distribution amplitude and the logarithmic corrections above the 
factorization scale correspond to the lower (upper) solid curves in the spacelike (timelike) 
regions. Similarly, the lower (upper) dashed curves correspond to our prediction in the 
spacelike (timelike) regions when using the asymptotic distribution amplitude. The param- 
eter /i was fixed from the known value of the normalization of the T form factor at Q 2 = 0. 
This is determined from the observed decay width it — > 27, T(0) ~ 0.14GeV _1 , which gives 
~ 338 MeV. This is certainly a mass scale typical of both confinement and dynamical 
chiral symmetry breaking as anticipated. The dotted curve is a simple dipole fit to the data 
using the vector meson dominance model (VMD), where T(Q 2 ) = T(0)/[1 + (Q 2 /M 2 )] with 
the parameter choice M ~ 750 MeV. Again, at very large \Q 2 \ solid and dashed curves 
converge separately. In contrast to the pion form factor case, however, this convergence 
seems to be much faster, and as shown Figs. |]a) and b) the relative 0(1/Q 4 ) correction to 
the pQCD predictions, 

ST{Q2) ^ T pQCD (Q 2 ) - T(Q 2 ) 
T p qcd{Q 2 ) 

is only about 30% (10%) and 20% (10%) for the Chernyak and Zhitnitsky distribution 
amplitude and asymptotic one in a spacelike (timelike) regions respectively. It is worth 
noting that the difference in the pQCD and our predictions for the T(Q 2 ) form factor below 
Q 2 < 50GeV 2 has opposite sign in the space- and timelike regions. This is due to the 
resonance structure in the timelike region which enhances the soft contributions and in 
particular introduces an imaginary part of the form factor. This feature is reproduced in 
our model and comes from the complex logarithm in Eq. (|37j) for Q 2 < —fi 2 . As — Q 2 
increases in the timelike region the two solid (dashed) curves cross and as —Q 2 — > 00 
they eventually merge as the imaginary part becomes negligible [i.e., ImT(Q 2 )/ReT(Q 2 ) — > 

log(-Q 2 /V)/ - Q 2 ] ■ 
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IV. CONCLUSIONS 



We have constructed a covariant, gauge- invariant qq amplitude for the pion in Minkowski 
space using a perturbative integral representation. We have shown that confinement can 
be implemented in a straightforward way. We have also shown that this wavefunction can 
be readily connected to the usual pQCD treatment in the asymptotic regime. The main 
conclusions of our pion electromagnetic form factor analysis come from an examination 
of the curves in Figs. |3] and [|. Comparing our predictions with the pQCD ones we see 
that the influence of the soft pion wavefunction on the hard gluon exchange diagrams can 
produce large effects. For example, our calculation of the soft contributions is as large 
as ~ 90% of the purge pQCD result for Q 2 ~ 5GeV 2 when the CZ distribution is used. 
For the asymptotic distribution our calculation of the soft contribution is ~ 50%. The 
incorporation of nonperturbative effects into a hard quark scattering calculation may be 
somewhat model dependent, but it seems clear that the effects can be very significant at 
intermediate momentum scales. Previous treatments of this subject ||, dealt only with 
the pure pQCD formula of Eq. ( p6|) and in particular they manifestly break the hard-soft 
factorization as the size of the soft corrections is estimated by modifying the hard scattering 
amplitude by a soft gluon mass. On the other hand, in our approach we study the sensitivity 
of the form factors to the soft physics coming from the wavefunctions while keeping the hard 
amplitude perturbative. In addition, the present study has been an attempt to partially 
address the issues of Lorentz covariance, gauge invariance, and the nonperturbative structure 
of the covariant amplitude of Eq. (|2|). Our results for the pion e.m. form factor appear to 
support the conclusion that the assumption of dominant pQCD contributions to the pion 
form factor at moderate Q 2 is not well-founded. 

We have also shown that the covariant pion wavefunction can be directly related to the 
matrix element associated with the 7*7r° — > 7 form factor. We have analyzed the extrapola- 
tion of this form factor away from the asymptotic region using the nonperturbative covariant 
pion wavefunction. Once we fit T(0), using the single nonperturbative scale parameter /1, 
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we find that T{Q 2 ) displays a much weaker dependence on the detailed form of the pion 
wavefunction even at low to moderate Q 2 , (i.e., \Q 2 \ < 10 GeV 2 ). In particular, we found 
that the relative, wavefunction-dependent, soft corrections to T(Q 2 ) are less than ~ 25% 
for \Q 2 \ ~ 10 GeV 2 and become even smaller as \Q 2 \ increases. Because of that and be- 
cause the predictions for the form factor for different types of distribution amplitudes lead 
to different predictions, a precision measurement of the forward angle cross section for the 
e+p — > e+p+7 transition would be of considerable interest, (e.g. at CEBAF). In particular, 
it may help resolve the question of large peaks predicted by Chernyak and Zhitnitsky for 
x -> ±1. 

The apparently conflicting conclusions drawn from the results for the pion e.m. form 
factor F(Q 2 ) and the 7* + tt° — > 7 from factor T(Q 2 ) require cautious interpretation. The 
lesson to be learned from these simple model analyses may be that leading pQCD treatments 
can be relatively successful in describing exclusive amplitudes in the intermediate momentum 
regime based on the operator product expansion [T(Q 2 )], but fail quite badly when there 
is no appropriate short distance expansion involved [e.g., F(Q 2 )} ||. Further studies along 
these lines will be needed before more definitive answers are possible. 
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FIGURES 



FIG. 1. An exact representation of the amplitude for the pion e.m. form factor. 

FIG. 2. Explicit one-gluon contributions to the high-Q 2 region of the form factor. The pQCD 
expression of Eq. (|l|) is obtained by multiplying the wave function by the collinear projectors, a) 
Hard gluon exchange at the soft pion vertices, b) spectator gluon (an example spectator diagram). 

FIG. 3. Spacelike region of the pion e.m .form factor. The various theoretical curves are 



explained in the text. The data are taken from Ref. [19] 



FIG. 4. Relative size of the nonperturbative contributions to the pion e.m. form factor for the 
CZ (solid) and asymptotic (dashed) distributions. 



FIG. 5. Amplitude for the j*tt — > 7 from factor 

FIG. 6. Diagrams contributing to the large-|Q 2 | behaviour of the T(Q 2 ) form factor. The 
pQCD prediction is obtained by multiplying the wavefunction by a collinear projection. 

FIG. 7. The j*ir — ► 7 transition form factor for a) spacelike and b) timelike virtual photon 
momenta. The various theoretical curves are explained in the text. The data are also shown |)23| . 

FIG. 8. An illustration of the relative magnitude of the soft contributions for the j*tt — > 7 
transition form factor in the a) spacelike and b) timelike regions. Again the solid and dashed 
curves are for the CZ and asymptotic distributions respectively. 
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